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Self-Diffusion  Coefficients  of  Sol-Gel  Intermediates 

L.  V.  Ng  and  A.  V.  McCormick*,  Chemical  Engineering  and  Materials  Science  Depanment, 
University  of  Minnesota,  Minneapolis,  MN  55455. 

ABSTRACT 

Over  the  past  few  decades,  the  sol-gel  technique  has  shown  great  potential  for 
producing  ceramic  materials  of  designed  composition  and  properties.  A  typical  sol-gel  process 
involves  the  hydrolysis  and  polymerization  of  one  or  more  metal  alkoxides  in  an  alcohol 
solution  in  a  batch  reactor  to  form  a  homogeneous  gel,  which  is  then  further  treated  to  provide 
the  type  of  ceramic  desired  (aerogel,  gi  iss,  crystal,  coating,  etc.). 

Several  authors  have  suggested  that  diffusion  effects  may  influence  the  development  of 
the  gel  microstructure.  To  investigate  if  diffusion  effects  arc  important,  we  begin  by 
determining  what  governs  the  diffusivity  of  the  precursors  using  the  pulscd-gradient-spin-echo 
NMR  technique.  This  report  focuses  on  the  early  stages  of  the  batch  reaction. 

The  system  chosen  is  the  acid  catalyzed  hydrolysis  and  condensation  of  tetraethyl 
orthosilicate  (TEOS)  in  ethanol.  Our  results  from  POSE  NMR  and  ^^Si  NMR  show  that 
the  diffusivities  of  the  silicate  precursors  are  not  just  a  function  of  the  molecular  weight 
Hydrogen  bonding  between  the  polar  functional  groups  and  the  solvent  also  decreases  the 
mobility  of  the  reactants.  Moreover,  since  hydrolysis  increases  and  polymerization  reduces 
the  concentration  of  these  polar  groups,  the  mobility  of  the  reactive  oligomers  changes  in  a 
complex  manner.  Combining  our  results  from  POSE  NMR  and  gas  chromatography,  we 
estimate  the  extent  of  hydrogen  bonding  effect  on  the  mobility  of  various  oligomers,  providing 
a  model  to  predict  diffusivity  trends  throughout  the  batch  reaction. 
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INTRODUCTION 

The  hydrolysis  and  condensation  of  TEOS  in  alcohol  is  one  of  the  most  extensively 
studied  systems  in  sol-gel  synthesis.  Previous  work  has  focused  on  the  reaction  mechanisms, 
the  kinetics,  the  structure  and  mechanism  of  growth,  and  the  effects  of  processing  parameters 
( for  an  excellent  review,  see  ref.  1 ).  Typical  reactions  are: 

Hydrolysis 

(EtO)3Si-OEt  H2O  (EtO)3Si-OH  +  EtOH  (1) 

Condensation 

(EtO)3Si-OEt  +  HO-Si(OEt)3  -»  (Et0)3Si-0-Si(0Et)3  +  EtOH  (2) 

(EtO)3Si-OH  +  HO-Si(OEt)3  ->  (Et0)3Si-0-Si(0Et)3  +  H2O  (3) 

where  Et  represents  CH2CH3.  The  reactions  above  represent  only  a  few  generic  reactions. 
The  entire  sol-gel  synthesis  process  involves  many  reactants  and  products  ranging  from 
monomers  to  large  polymers,  each  with  various  degrees  of  hydrolysis. 

X-ray  scattering  studies  [  2-5  ]  have  shown  that,  depending  on  processing  conditions 
( e.g.  type  of  catalyst,  water  concentration,  etc. ),  the  condensation  reactions  can  lead  to 
radically  different  gel  structures  -  from  colloidal  particles  to  weakly  branched  polymeric 
networks.  Results  obtained  by  Himmel  and  coworkers  [  3  ]  show  that  the  fractal  dimension 
of  the  gels  decreases  with  decreasing  starting  water  concentration  and  solution  pH  ( for  pH 
values  less  than  2  ).  Consistent  with  other  results  [  4,5  ],  gels  produced  with  low  water 
concentration  at  low  pH  were  found  to  be  weakly  branched. 

Computer  simulations  [  6  ]  on  the  growth  of  fractal  objects  have  shown  that  the  fractal 
dimension  can  depend  on  the  nature  of  the  growth  process  ( e.g.  particle-cluster  vs.  cluster- 
cluster  growth )  and  the  rate-determining  step  (  e.g.  reaction-limited  vs.  diffusion-limited ). 
Comparing  their  results  with  con^uter  simulations,  Himmel  and  others  suggest  that  the  weakly 
branched  gel  produced  from  the  hydrolysis  and  condensation  of  TEOS  with  low  water 
concentration  and  solution  pH  is  a  result  of  a  diffusion-limited  growth  process. 
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Pouxviel  and  Boilot  [  7  ]  have  found  that  even  though  initially  the  system  with  water- 
to-silicon  molar  ratio  ( )  of  10  showed  faster  hydrolysis  rate  and  higher  degree  of 
condensation  than  the  system  with  =  4,  the  former  takes  about  twice  as  long  to  gel.  To 
explain  this  apparent  inconsistency  in  kinetic  behavior,  they  speculate  that  the  growth  process 
changes  from  kinetic-  to  diffusion-controlled  as  the  polymers  become  sufficiently  large,  at 
which  point  the  diffusion-controlled  growth  is  then  faster  for  the  system  with  smaller  polymers 
( i.e.  Fvv  =  4 ).  Colby  and  coworkers  [  8  ]  compared  the  gel  rime  of  an  HCl-catalyzed 
TEOS  system  with  a  similar  TMOS  ( tetramethyl  orthosilicate  )  system,  and  found  that  the 
apparent  activation  energy  for  gelation  was  higher  for  the  TEOS  system.  The  difference  was 
attributed  by  these  workers  to  the  effects  of  diffusion. 

The  origin  of  different  gel  structures  is  still  not  well  understood.  In  panicular,  it  is  not 
clear  how  diffusion  limitations  can  play  a  role  in  the  microstructural  development  of  the  gels. 
By  understanding  factors  controlling  diffusion,  it  is  possible  to  engineer  the  properties  of  sol- 
gel  derived  materials.  In  this  work,  we  investigate  the  factors  that  govern  the  mobility  of 
intermediates.  It  is  to  be  expected  that  changes  in  diffiisivity  might  be  more  complicated  than 
in  ordinary  polymer  solutions  since  molecular  weight,  solvent  polarity,  and  polymer  group 
polarity  are  changing  simultaneously  in  a  batch  reactor. 
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EXPERIMENTS 

1 .  Preparation  of  the  TEOS  reaction  mixture 

TEOS  and  absolute  ethanol  (  both  used  as  purchased  from  Aldrich  and  Minnesota 
Grain  Products,  respectively )  were  mixed  slowly  for  20  seconds.  Aqueous  acid  solution 
was  added,  with  vigorous  stirring,  for  1  minute  to  yield  a  solution  with  the  molar  composition 
TEOS  :  EtOH  :  H2O  :  HCl  =  1:4:  0.8  :  0.0008  ( the  selection  of  this  system  will  be 
discussed  later).  The  sample  was  kept  at  room  temperature  and  remained  clear  over  the 
reaction  times  studied  ( from  1  hour  to  27  days ).  Samples  were  prepared  with  a  range  of 
reaction  times  and  were  analyzed  together  on  the  same  day  so  that  signal  intensities  could  be 
compared  quantitatively  between  spectra  without  variations:  introduced  by  probe  tuning,  magnet 
shimming,  etc..  Three  separate  sets  of  data  were  obtained  to  estimate  the  reproducibility  of  thf* 
diffusion  results. 

2.  iH  and  29si  nMR 

spectra  were  obtained  on  a  Nicolet  NT300  at  300  MHz  and  spectra  on  a 
Varian  VXR500  at  99.3  MHz.  Chromium(III)  acetylacetonate  was  used  as  a  paramagnetic 
relaxing  agent  to  reduce  the  Tj  of  29si  to  about  2.5  seconds.  For  each  ^^Si  spectrum,  a 
total  of  1024  scans  separated  by  a  relaxation  delay  of  3  seconds,  and  a  line- broadening  of 
1.0  Hz  were  employed.  To  check  if  the  fast  relaxation  delay  affect  the  quantitative  nature  of 
the  spectra,  we  compared  the  integrated  peak  intensities  of  these  spectra  with  those  where  an 
ample  relaxation  delay  of  12  seconds  was  used  ( Figures  8  and  13  );  we  have  found  no 
significant  difference. 

3.  PGSE  NMR  experiment 

The  pulsed-gradient-spin-echo  NMR  technique  is  used  to  determine  the  self-diffusion 
coefficient  of  inteimediates  oligomers.  The  technique  involves  labeling  the  positions  of  nuclei 
by  applying  a  magnetic  rield  with  a  spatial  gradient  It  is  a  non-intrusive  method  and  ideal  for 
probing  the  thermal  motion  of  chemical  species.  The  technique  is  standard  for  studies  in 
polymer  solutions  and  melts  [e.g.,  9-12]  and  emulsions  [e.g.,  13].  Detailed 
descriptions  of  the  PGSE  NMR  technique  can  be  found  in  several  sources  [  9,10,13,14  ], 
which  explain  that  the  self-diffusion  coefficient  D,  of  a  chemical  species  is  related  to  its  NMR 
signal  intensity,  S,  by  the  following  equation  [  10  ]: 


5  =  Soexp(-2A/T2)exp(-72G2|3D)  (4) 

where. 

|3  =  52(  A  -  6/3  ),  (5) 

T2  =  the  spin-spin  relaxation  time  constant 
G  =  the  field  gradient  strength,  in  units  of  Tesla»cm'* 

A  =  the  diffusion  observation  time 

6  =  the  duration  of  the  pulsed  gradient 

Y  =  the  nucleus  gyromagnetic  ratio,  in  units  of  rad*Tesla‘  1,5-1 

PGSE  experiments  were  performed  on  a  Nicolet  300  MHz  NMR  spectrometer  with 
parameters: 

A  =  100  ms 

6  =  0.4  ms  to  5.4  ms 

G  -  15  Gauss/cm  =  0.15  Tcsla/m 

For  self-diffusion  coefficients  in  the  order  of  10'^  cm2/s,  an  observation  time  of  A  =  100  ms 
is  recommended  [  15  ].  This  corresponds  to  a  root-mean-square  displacement  of: 


The  field  gradient  strength  was  calibrated  using  benzene  at  25  °C  ( 2.28  x  10'^  cm2/s ). 
For  each  self-diffusion  measurement,  10  to  26  values  of  6  were  used,  holding  A  constant 
In  accordance  with  equation  4,  the  self-diffusion  coefficient  D  of  the  silicate  oligomers  was 
found  by  plotting  ln(S)  versus  P,  where  S  is  the  integrated  proton  signal  intensity  of  the 
(-SiO-CH2)  group. 

4.  Gas  Chromatography 

An  experimental  study  of  imlecular  motion  in  a  multi-component  system  must  include 
the  identification  of  various  molecules  and  the  determination  of  their  concentrations.  These 
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molecules  are  the  various  silicate  oligomers  ( e.g.  monomers,  dimers,  trimers,  etc. ),  the 
distribution  of  which  cannot  be  accurately  determined  solely  by  NMR.  Combining  data  from 
gas  chromatography  and  ^^Si  NMR,  though,  we  were  able  to  separate  and  identify  the  silicate 
oligomers.  These  techniques  have  been  used  by  other  investigators  [  16-18  J  for  studies  on 
the  hydrolysis  and  condensation  of  TMOS  ( tetramethyl  orthosilicate  ). 

Flame  ionization  detection  was  effective  in  detecting  and  quantifying  hydrocarbon- 
containing  molecules  [  19  ].  It  has  been  referred  to  as  a  carbon-counting  device  since  the 
response  factor  for  a  molecule  is  approximately  proportional  to  the  number  of  carbon  it 
contains  ( comparison  of  GC-FID  peak  intensities  with  ^^Si  NMR  results  showed  that  the 
approximation  is  reasonable ).  We  were  able  to  quantify  the  molecular  size  distribution  using 
this  technique  since  our  system  does  not  contain  completely  hydrolyzed  silicates,  as  verified  by 
29si  NMR. 

The  GC-FID  used  was  a  Hewlett  Packard  HP  5890  series  II.  A  10  m  x  0.53  mm 
i.d.  fused  silica  capillary  column  bonded  with  50%  phenyl  and  50%  methyl-polysiloxane,  film 
thickness  2.0  |im  ( Hewlett  Packard  model  HP- 17  ),  is  employed  using  helium  as  carrier  gas. 
The  operating  conditions,  following  [  17  ],  were: 

injector  temperature  =  300 
detector  temperature  =  300  <^3 
oven  temperature  program: 

60  oC  for  1  minute 
15  0C/min  to  150  oC 
5  0C/min  to  260  «€ 

260  °C  for  1  minute, 
column  flow  rate  =  22  ml/min 


Figure  4 
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RESULTS  AND  DISCUSSION 


A.  Selection  of  Composition 

A  typical  proton  NMR  spectrum  for  the  reaction  mixture  is  shown  in  Figure  1.  Note 
that  because  of  fast  chemical  exchange,  the  (-OH)  peak  is  a  mole-weighted  average  of  the 
signals  from  the  protons  of  water,  the  hydroxyl  groups  on  ethanol  and  silanols.  The 
phenomenon  of  fast  chemical  exchange  occurs  because  the  nuclei  are  rapidly  exchanging 
between  various  molecular  sites,  with  respect  to  the  NMR  detection  time.  As  the  exchange 
slows  down,  the  individual  peaks  will  be  resolved.  This  effect  can  be  seen  on  a  spectrum 
of  a  mixture  of  TEOS,  ethanol  and  water  ( Figure  2 ).  The  exchange  slows  down  without  the 
acid,  which  catalyzes  proton  exchange  [  20  ].  Note  that  a  spectrum  (  Figure  3  )  with 
higher  water  concentradon  confirms  the  water  peak  idendficadon.  High  water  concentradon 
promotes  proton  exchanges,  as  shown  by  the  collapse  of  ethanol  hydroxyl  triplet. 

Another  feature  of  the  spectra  is  that  protons  on  all  silicon-bonded  ethoxy  groups 
resonate  at  virtually  the  same  frequency.  This  means  that  various  silicate  species,  e.g. 
monomers,  dimers,  etc.,  cannot  be  differendated  through  the  proton  signals  on  the  NT300  ( a 
300  MHz  spectrometer )  [21  ].  Moreover,  the  methyl  signal  on  the  silicate  overlaps  with 
that  on  the  ethanol.  Since  the  (•SiO-CH2)  signal  is  the  only  one  available  for  the  labeling  of 
the  various  silicate  species,  the  self-diffusion  coefficient  measured  is  an  average  of  all  ethoxy- 
containing  silicates.  In  order  to  preserve  this  signal  for  the  PGSE  experiments,  we  are 
restricted  to  systems  with  low  concentradon  of  water  because,  when  present  in  high 
concentration,  water  promotes  hydrolysis  which  replaces  the  ethoxy  groups,  as  shown  in 
Figure  4.  These  systems  also  tend  to  yield  weakly-branched  polymers  which,  simuladons 
suggest,  are  formed  by  a  diffusion-limited  process.  Further  work  is  underway  using 
"tagging"  silyladon  technique. 
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B .  Self'diffusion  of  silicate  intermediates 

F’.gure  5  shows  that  the  plot  of  ln(S)  versus  P,  from  which  the  average  self-diffusion 
coefficient  of  the  ethoxy-containing  silicates  is  obtained,  is  a  straight  line,  suggesting  that 
polydispersity  in  the  silicate  mobility  is  negligible  [  22  ].  The  change  in  the  silicate  mobility 
with  reaction  is  shown  in  Figure  6. 

The  observed  initial  decrease  in  mobility  is  consistent  with  the  expected  increasing 
average  molecular  weight  of  the  silicates.  The  most  striking  feature  of  the  diffusion  data  is 
that,  rather  than  continuing  to  decrease,  the  average  mobility  of  the  ethoxy-containing  silicates 
levels  off,  and  then  rises  slightly,  even  though  we  expect  that  the  molecular  weight  continues  to 
increase.  This  anomaly  suggests  that  there  are  factors  other  than  the  molecular  weight  which 
determine  the  silicates  mobility. 

To  reveal  other  factors  affecting  the  mobility  of  the  silicates,  ^^Si  NMR  is  used. 


C.  The  effect  of  molecular  weight 
1.  29si  NMR 

The  peak  assignments  for  the  ^^Si  NMR  spectra,  displayed  in  Figure  7,  follow  those 
reported  in  the  literature  [  21(b),23,24  ].  The  chemical  shifts  of  commercial  products  TEOS 
and  ethyl  polysilicates  are  known.  The  latter  exhibits  three  lines  corresponding  to  the  end, 
middle  and  trisubstituted  groups.  In  addition,  since  a  hydroxyl  substituent  on  the  silicon 
nucleus  is  more  electron-withdrawing  than  an  ethoxy  substituent,  a  shift  of  -i-2.S  ppm  for  each 
replacement  of  an  ethoxy  group  by  a  hydroxyl  group  has  been  shown  [  21(b)  ].  The 
assignments  of  the  unhydrolyzed  middle-group  silicon  in  a  cyclic  tetramer  ( approx. 
-95  ppm )  and  a  linear  oligomer  ( approx.  -96  ppm )  follow  the  results  obtained  by 
Marsmann  and  coworkers  [  25  ].  The  strained  bond-angle  in  a  cyclic  tetramer  causes  a  slight 
deshielding  effect.  This  trend  is  also  observed  by  Harris  and  Knight  who  identified  the  ^Si 
chemical  shifts  of  the  silicate  ions  in  alkaline  aqueous  potassium  silicate  solutions  [  26  ]. 
Note  that  no  cyclic  trimers  are  found  ( the  unhydrolyzed  and  completely  hydrolyzed  species 
resonate  at  -87.97  ppm  [  25  ]  and  -82.9  ppm  (  27  ],  respectively  ). 
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The  multiplets  for  (SiO)-Si-(OH)(OEt)2,  (SiO)-Si-(OEt)3  and  linear  (SiO)2-Si-(OEt)2 
sites  are  due  to  the  variety  of  different  environment  of  the  adjacent  silicon.  Hydrolysis  of  the 
adjacent  silicon  causes  deshieldmg,  while  substitution  on  this  silicon  causes  shielding.  Peak 
intensity  comparisons  of  several  spectra,  combined  with  results  from  GC-FID,  lead  to  the  more 
detailed  assignments  shown  in  Figure  8.  These  assignments  agree  with  those  reported  by 
Kelts  and  Armstrong  [  24  ]. 

It  can  be  seen  from  the  ^^Si  NMR  spectra  (  Figure  7  )  that  all  the  silicons  are  linked  to 
ethoxy  groups,  so  that  they  are  accounted  for  in  the  diffusion  and  GC-FID  measurements.  To 
examine  the  evolution  of  different  silicon  sites,  the  peaks  are  integrated,  and  the  intensities 
plotted  against  reaction  time  in  Figure  9. 

From  Figure  9,  it  appears  that,  between  0  and  20  hours,  a  high  concentration  of  the 
monohydrolyzed  monomers  ( Si(OEt)3(OH) )  is  formed,  and  quickly  condensed,  resulting  in 
high  concentration  of  endgroups  (SiO)-Si(OEt)3.  A  combination  of  hydrolysis  and 
condensation  also  results  in  (SiO)-Si(OEt)2(OH).  These  hydrolyzed  endgroups  undergo 
further  condensation,  leading  to  the  disubstituted  silicons  ( (SiO)2-Si(OEt)2  ).  The  liberated 
water  hydrolyzes  the  remaining  endgroups,  resulting  in  further  condensation.  The  reactions 


can  be  represented  as  ( R  can  be  Et  or  H ): 

Si(OEt)4  +  H2O  (HO)-Si  (OEt)3  +  EtOH  (7) 

sSi-OR  +  HO-Si-(OEt)3  ->  sSi-0-Si(0Et)3  +  ROH  (8) 

sSi-OR  +  (HO)2-Si-(OEt)2  ->  HSi-0-Si(0H)(0Et)2  +  ROH  (9) 

sSi-0-Si(0H)(0Et)2  +  RO-Sis  ^  (HSi-0)2-Si(0Et)2  +  ROH  (10) 


When  the  concentration  of  water  and  silanols  drop  ( time  >  200  hours  ),  the  rates  of 
hydrolysis  and  condensation  should  become  comparable.  Consequently,  as  shown  by  the 
following  equations,  the  concentration  of  silanol  decreases  only  slowly. 


sSi-OH  +  HO-Sia  sSi-O-Sis  +  H2O 
aSi-OEt  +  H2O  ->  sSi-OH  +  EtOH 


(11) 

(12) 
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It  is  clear  from  Figure  9  that  condensation  continues  throughout  the  experiments,  as 
shown  by  the  increase  in  concentration  of  the  disubstituted  silicon  groups  and  by  the 
decreasing  proportion  of  (SiO)-Si(OEt)2(OH).  Since  the  (SiO)-Si(OEt)3  sites  are  present  in 
much  higher  concentration  than  the  (SiO)-Si(OEt)2(OH)  sites,  the  condensation  observed  here 
will  most  likely  result  in  short  chain  oligomers  ( chances  of  termination  is  high  ).  Moreover, 
since  the  degree  of  hydrolysis  is  low,  the  oligomers  are  probably  linear.  These  observations 
are  supported  by  the  results  of  GC-FID. 

2 .  Gas  Chromatography 

Gas  chromatograms  are  shown  in  Figure  10,  where  the  peak  assignments  are  a^ed  by 
the  results  of  Klemperer  and  coworkers  [  17,28  ]  as  well  as  the  “^Si  NMR  results.  Similar 
oligomers  have  been  detected  for  systems  with  the  same  range  of  water-to-silicon  ratio  [  29  ]. 
It  has  been  found  [  28  ]  that  the  presence  of  hydroxyl  groups  on  an  oligomer  increases  the 
species  retention  time  slightly.  This  is  understandable  since  the  hydroxyl  groups  promote 
hydrogen  bonding  which  will  lower  the  vapor  pressure. 

The  mole  fractions  of  the  silicon  oligomers  are  calculated  from  the  peak  intensities,  and 
are  shown  in  Figure  11.  These  results  show  that  reactions  are  occurring  rapidly  within  the 
first  few  hours  after  mixing,  leading  to  the  sequential  formation  of  hydrolyzed  monomers, 
dimers  and  trimers.  The  reaction  rates  slow  markedly  after  24  hours,  but  condensation  is  still 
occurring.  Note  that  beyond  a  certain  point,  all  reactions  have  ceased,  so  that  there  is  no 
observable  change  between  the  gas  chromatograms  for  the  4-month  and  9-month  samples. 


D.  Hydrogen  Bonding 
1 .  Small  oligomers 

It  is  clear  that  some  factor  other  than  molecular  weight  is  decreasing  the  average 
mobility  of  the  silicates.  The  solvent  is  a  mixture  of  water  and  ethanol,  containing  polar 
hydrogen  atoms  which  are  ready  to  form  hydrogen  bonds  with  the  hydroxyl  groups  of  the 
silanols.  The  mobility  of  these  solutes  will  consequently  decrease  [  30-32  ].  The  effect  of 
hydrogen  bonding  on  the  solute  mobility  has  been  reported  for  some  polymeric  systems 
[  33  ]. 
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One  indication  of  the  H-bonding  effect  is  the  drastic  drop  in  the  silicate  mobility  one 
hour  after  mixing,  when  chromatograms  show  that  low  molecular  weight  silanols  are  present  in 
high  concentrations.  After  24  hours,  GC-FID  and  ^^Si  NMR  (  Figure  9  )  both  show  that 
the  reaction  mixture  still  contains  a  high  concentration  of  silanols,  so  the  silicate  mobility  stays 
rather  constant  owing  to  the  solute-solvent  associations  effect.  Eventually,  when  the  silanols 
are  consumed  completely,  a  rise  in  silicate  mobility  is  observed.  The  loss  in  silicon  hydroxyl 
groups  is  also  coincidental  with  the  complete  consumption  of  water  (  Figure  12  ).  Without 
any  water,  reactions  stop  at  this  point,  leaving  predominantly  dimers  and  trimers 
( Figure  10  )  with  an  average  diffusion  coefficient  of  0.60  x  lO"^  cm^/s. 

The  consumption  of  water  is  a  net  result  of  the  hydrolysis  and  condensation  of  TEOS, 
as  shown  by  the  following  stoichiometric  relations: 

2  (EtO)3Si-OEt  +  2H2O  ^  2(EtO)3Si-OH  -t-  2EtOH  (13a) 

2  (EtO)3Si-OH  ^  (Et0)3Si-0-Si(0Et)3  +  H2O  (13b) 

The  equation  above  indicates  that  water  is  gradually  consumed  and  ethanol  produced  as 
reactions  progress.  The  system  starts  with  low  water  concentration  ( a  substoichiometric 
water-to-silicon  molar  ratio,  W/Si,  of  0.8  ),  so  that  the  water  supply  for  the  reactions  vanishes 
before  an  appreciable  degree  of  polymerization  is  achieved.  It  is  shown  in  Appendix  A  how 
the  water  and  ethanol  content  of  the  system  can  be  calculated  from  the  ‘^Si  NMR  results  using 
molar  balances. 

The  effect  of  H-bonding  in  reducing  the  mobility  of  silicon  precursors  can  no  doubt  i  ^ 
better  demonstrated  by  studying  systems  with  high  degree  of  hydrolysis.  However,  this  plat; 
is  limited,  at  present,  by  experimental  considerations  since  we  must  retain  the  silicon-bonded 
ethoxy  group  for  NMR  labeling.  Even  for  systems  which  undergo  limited  hydrolysis,  the  rate 
of  reaction  must  be  slow  enough  so  that  diffusion  measurements,  which  typically  runs  for 
about  7  minutes  for  one  result,  can  be  performed.  However,  there  is  hope  of  "tagging"  large, 
fully  hydrolyzed  polymers  with  a  trimethylsilylating  reagent  [  34  ]. 

To  examine  the  effect  of  H-bonding  more  closely,  the  extent  of  hydrolysis  on 
intermediate  silicon  oligomers  can  be  increased  by  first  preparing  incompletely  hydrolyzed 
oligomers  and  then  hydrolyzing  them  with  additional  water. 
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2.  Intermediate  oligomers 

As  shown  earlier,  a  system  with  low  water  concentration  produces  only  silicon 
oligomers,  all  containing  ethoxy  groups.  To  such  a  system,  addition  of  limited  amount  of 
water  at  the  end  of  the  reaction  should  again  initiate  hydrolysis  but  without  losing  all  the  ethoxy 
groups  for  diffusion  measurements.  Moreover,  since  hydrolysis  is  occurring  with  silicon 
oligomers,  the  rate  of  reaction  is  slower  than  if  it  is  with  monomers  [  7,35  ],  and  the 
conditions  are  more  closely  related  to  those  in  our  previous  experiment 


To  produce  the  initial  silicon  oligomers,  a  system  with  molar  composition  TEOS  : 
EtOH  :  H2O  :  HCl  =  1  :  4  :  0.5  :  0.0005  is  prepared.  The  sample  is  then  allowed  to  react  for 
about  3  months  ( the  29$  i  NMR  spectrum,  shown  in  Figure  13,  for  a  1 -month-old  sample 
shows  that  reactions  have  virtually  ceased,  leaving  monomers,  dimers  and  trimers ).  To  the 
3-month-old  sample,  limited  amounts  of  water  ( total  W/Si  =  0.8  and  1.0  )  are  added,  and 
ordinary  NMRandPGSE  NMR  experiments  are  performed.  To  estimate  the  extent  of 
hydrolysis  (  without  resorting  to  29si  NMR ),  the  molar  ratio  of  the  silicon-bonded  ethoxy 
f  -SiO-CHo  -  1 

groups  to  ethanol  hO-Ch/  J  **  NMR  spectra.  Hydrolysis 

leads  to  a  drop  in  the  concentration  of  the  ethoxy  group  and  a  rise  in  ethanol  concentration. 

”~SiO-“CH*—  "1 

Consequently,  the  extent  of  initial  hydrolysis  is  reflected  by  the  drop  in  the  | 

HO 


- - - -  - . - .  -  - r - 1^  HO-CH2-  J 

molar  ratio.  This  ratio  levels  off  when  hydrolysis  slows  down  and  condensation  becomes 
significant,  and  so  this  ratio  characterizes  the  amount  of  hydrolyzed  groups  present  Hence, 

r  -SiO-CH?-  1 

by  monitoring  and  comparing  the  ---jr  molar  ratio  and  the  silicate  average 

L  Hi>-CH2—  j 

mobility,  the  effect  of  H-bonding  can  be  examined  ( Figure  14 ). 


molar  ratio  and  the  silicate  average 


The  initial  drop  in  the  average  mobility  follows  the  extent  of  hydrolysis  consistently  for 
both  systems.  The  initial  larger  drop  in  mobility  for  the  system  with  higher  water 
concentration  is  consistent  with  a  higher  degree  of  hydrolysis.  As  hydrolysis  slows  down  and 
condensation  speeds  up,  the  average  mobility  rises  -  a  trend  similar  to  the  one  observed  in 
previous  experiment.  Notice  that  with  the  loss  of  hydrogen  bonding,  the  difference  in 
mobility  between  the  two  systems  becomes  smaller.  This  indicates  that,  at  least  for  low 
molecular  weight  systems,  the  effect  of  hydrogen  bonding  is  more  significant  than  the 
molecular  weight  effect 
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Results  from  the  diffusion  experiments  have  shown  that  the  diffusivity  of  silicon 
precursors  is  reduced  both  by  the  increase  in  their  molecular  weight  and  the  degree  of 
hydrolysis,  the  latter  is  due  to  hydrogen  bonding  with  the  solvent  molecules.  The  extent  of 
these  two  factors  can  be  seen  by  comparing  the  estimated  diffusivities  of  the  hydrolyzed  and 
unhydrolyzed  species. 


E.  Diffusivity  of  Hydrolyzed  and  Non-Hydrolyzed  Oligomers 

The  average  diffusion  coefficient  measured  by  *H  POSE  NMR  experiments  at  any 
time  is  related  to  the  self-diffusion  coefficient  of  the  silicon  oligomers  by  the  following 
equation: 


m 

^av  - 

j=l 

cth 

where  fj  is  the  fraction  of  ethoxy  belonging  to  the  j-mer  ( the  oligomer  containing  j 
silicon )  and  can  be  obtained  from  the  gas  chromatograms  since  the  GC-FID  peak  intensity  is 
proportional  to  the  number  of  carbon  in  the  molecule.  Thus, 

eth  (%Area)i 
■  m 

^(%Area)ij 

k»l 


where  (%Area)j  is  the  percent  integrated  peak  area  for  species  j.  Between  10  and  40 
hours,  the  species  that  are  present  in  signifrcant  amount  are  monomer,  monohydrolyzed 
monomer,  dimer,  monohydrolyzed  dimer,  trimer  and  monohydrolyzed  trimer.  To  extract  the 
diffusion  coefficients  of  these  species  from  the  average  value,  we  need  to  solve  a  set  of  6-by-6 
linear  equations.  Since  we  presume  that  the  monomer  diffusivity  is  the  average  diffusivity 
obtained  at  time  t  =  0,  we  only  need  S  more  data  sets.  The  percent  peak  area  in  equation  (IS) 
and  the  average  diffusivities  are  extracted  from  the  curves  drawn  through  the  corresponding 
data.  The  linear  equation  set  is  solved  using  the  method  of  Gaussian  elimination  with  partial 
pivoting  [  36  ],  and  the  results  are  shown  in  Table  1  and  plotted  in  Figure  IS.  The  Stokes- 
Einstein  hydrodynamic  radii  are  estimated  using  a  pure  ethanol  solvent  viscosity  of  1.1  cP 
obtained  from  the  literature  [  32  ]. 


Oligomer 


Diffusivity  ^  [  10'^  cirr/s  ] 


monomer 

0.900 

2.206 

dimer 

0.591 

3.360 

trimer 

0.521 

3.810 

linear  tetramer  ^ 

0.500 

3.971 

monohydrolyzed  monomer 

0.447 

4.441 

monohydrolyzed  dimer 

0.176 

11.283 

monohydrolyzed  trimer^ 

0.100 

19.850 

^  The  error  is  estimated  to  be  ±  0.03  based  on  the  deviation  on  the  average  diffusivity. 
^  From  extrapolation. 


Note  that  the  diffusivities  above  are  estimated  based  on  the  average  diffusivities  measured 
when  the  water  concentration  is  still  high.  Complexation  with  the  water  molecules  reduces  the 
diffusivities  of  the  hydrolyzed  species  by  increasing  the  effective  hydrodynamic  radii.  At  this 
stage  of  the  reaction,  molecular  weight  plays  a  lesser  role  in  determining  the  diffusivity  of  the 
oligomers. 

The  hydrogen  bonding  effect  in  reducing  the  mobility  of  the  silicon  species  has  been 
demonstrated.  Nevertheless,  the  consequences  of  H-bonding  are  far-reaching.  The  effects  of 
H-bonding  on  chemical  reaction  rates  have  been  observed  [  37  ].  Moreover,  according  to  the 
absolute  reaction  rate  theory  [  38  ],  H-bonding  can  reduce  the  condensation  reaction  rate 
constant  by  increasing  the  heat  of  activation.  This  effect  suggests  that  condensation  is  more 
favorable  in  aprotic  solvents.  It  also  suggests  that  alcohol-producing  condensation  reaction  is 
thermodynamically  more  favorable  than  water-producing  condensation  since  the  latter  involves 
two  water-complexed  reactants.  The  effects  of  hydrogen  bonding  on  the  condensation  process 
in  sol-gel  reactions  have  been  discussed  in  the  literature  [  39  ]. 

Comparing  our  diffusion  results  with  the  kinetic  studies  performed  separately  in  our 
group,  the  Smoluchowski  equation  [  40  ]  predicts  that  the  early  stages  of  condensation 
reactions  ( i.e.  when  the  Brownian  motion  of  the  reactants  is  not  hindered )  are  reaction- 
limited.  However,  when  the  polymers  become  sufficiently  large,  diffusion  effects  can  be 


important.  The  influence  of  diffusion  on  the  overall  reaction  fate  and  the  polymer 
microstructural  development  is  presently  being  examined 


SUMMARY 


It  has  been  suggested  in  the  sol-gel  literature  that  diffusion  effects  may  influence  the 
development  of  the  gel  microstructure.  To  investigate  if  diffusion  effects  are  important,  we 
begin  by  determining  what  governs  the  diffusivity  of  the  precursors 

This  work  focuses  on  the  early  stages  of  the  sol-gel  reaction  by  studying  the  acid- 
catalyzed  hydrolysis  and  condensation  of  TEOS  with  starting  molar  ratio  TEOS  :  Ethanol 
;  H2O  :  HCl  =  1  :  4  :  0.8  :  0.0008  using  ^H,  29si  and  ‘H-PGSE  NMR,  and  gas 
chromatography  with  flame  ionization  detector.  Low  molecular  weight,  partially  hydrolyzed 
oligomers  ( mainly  dimers )  are  formed  rapidly  and  the  average  mobility  of  these  silicates 
decreases  from  0.90  x  10*^  cm^/s  ( monomer  diffusivity )  to  an  average  value  of 
0.56  X  10'^  cm^/s.  As  condensation  continues,  shown  by  ^^Si  NMR  and  GC-FID,  the 
average  mobility  remains  around  this  value  up  to  250  hours  and  then  goes  up  slightly  to 
0.60  X  10'^  cm^/s.  Our  results  show  that  the  mobility  of  the  silicon  precursors  is  reduced 
with  increasing  molecular  weight  and  the  extent  of  hydrolysis.  Complexation  with  the  water 
molecules  reduces  the  diffusivities  of  the  hydrolyzed  species  by  increasing  the  effective 
hydrodynamic  radii.  At  early  stage  of  the  reaction,  molecular  weight  plays  a  lesser  role  in 
determining  the  diffusivity  of  the  oligomers.  We  attribute  the  late  rise  in  mobility  to  the  loss  in 
hydrogen  bonding  as  silanols  and  water  are  consumed  completely.  Reactions  have  stopped 
without  water  so  that  our  system  produces  only  low  molecular  weight  oligomers,  the  largest 
detected  in  significant  amount  being  the  linear  tetramers.  Although  the  early  stages  of 
reactions  are  clearly  not  diffusion-limited,  our  results  suggest  that  hydrogen  bonding  may  be 
important  in  slowing  the  motion  of  intermediates  in  the  presence  of  high  water  concentration. 

These  results  have  demonstrated  that  unlike  ordinary  polymer  solutions,  the  diffusivity 
of  sol-gel  intermediates  depends  on  a  combination  of  molecular  weight,  .solvent  polarity,  and 
polymer  group  polarity,  the  contribution  of  each  of  these  factors  changes  with  the  extent  of 
reaction. 
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Appendix  A 


Water  and  Ethanol  Content 

From  equations  (1),  (2)  and  (3),  we  observe  that,  in  the  hydrolysis  and  condensation  of 
TEOS,  the  proton  (H)  and  ethyl  (Et)  groups  are  neither  created  nor  destroyed.  Hence,  the  total 
concentration  of  each  of  these  groups  is  conserved.  A  molar  balance  on  (Et-)  shows  that  the 
ethanol  concentration  at  any  time  ( expressed  in  molar  ratio  of  ethanol-to-silicon )  is  related  to 
the  concentrations  of  various  silicon-groups  as  determined  from  the  ^^Si  NMR  spectra  in 
Figure  9  ( all  quantities,  in  bold  letters,  are  expressed  as  ratios  to  the  total  concentration  of 
silicons,  subscript  zero  indicates  initial  value  ): 

Et-OH  +  EtO-Si  =  (Et-)o  (A-1) 

EtO-Si  =  3  [(EtO)3-Si-OH]  +  4  [(EtO)4-Si] 

-h  2  [(EtO)2(OH)-Si-(OSi)]  +  3  [(EtO)3-SI-(OSi)]  +  2  [(EtO)2-Si-(OSi)2] 

(A-2) 

From  the  starting  composition,  the  initial  concentration  of  the  ethyl  group  is  8.  Similarly,  the 
water  concentration  at  any  time  can  be  determined  from  a  molar  balance  of  (H-),  as  shown 
below: 


H2O  +  EtOH  +  Si-OH  =  (H-)o  (A-3) 

Si-OH  =  [(EtO)3-Si-OH]  +  [(EtO)2(OH)-Si-(OSi)]  (A-4) 

The  value  of  EtOH  in  equation  (A-3)  can  be  obtained  from  equations  (A-1)  and  (A-2),  while 
the  initial  concentration  of  (H-)  is  5.6.  The  concentration  of  water  is  then  equal  to  j  [H2O]. 

Figure  12  shows  the  calculated  results. 
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FIGURE  CAPTIONS 


1.  A  typical  *H  spectrum  ( ref.  TMS  )  for  a  reaction  mixture  of  TEOS,  ethanol,  water 
and  HCl. 

molar  ratio  TEOS  :  EtOH  :  H2O :  HCl  =  1 : 4 : 0.8  : 0.0008 
reaction  time  =  655  hours  ( *27  days ) 

2.  A  *H  spectrum  for  a  mixture  of  TEOS,  ethanol  and  water, 
molar  ratio  TEOS  :  EtOH :  H2O  =  1 : 4 : 0.8 

The  spectmm  was  obtained  within  about  0.5  hour. 

No  reaction  could  have  occurred  appreciably. 

3.  A  ^H  spectrum  for  a  mixture  of  TEOS,  ethanol  and  water, 
molar  ratio  TEOS  :  EtOH :  H2O  =1:4:4 

The  spectrum  was  obtained  within  about  1  hour. 

No  reaction  could  have  occurred  appreciably. 

4.  A  ^H  spectrum  for  a  reaction  mixture  of  TEOS,  ethanol,  high  water  concentration  and 
HCl 

molar  ratio  TEOS  :  EtOH :  H2O  :  HC1=  1 : 4 : 4 : 0.0005 
reaction  time  =  26  hours 

High  water  concentration  results  in  significant  extent  of  hydrolysis,  and  thus  the 
disappearance  of  silicon-bonded  ethoxy  signals. 


5.  Determination  of  self-diffusion  coefficient  from  the  slope  of  a  plot  of  ln(S)  versus  p, 

where  S  is  the  (-SiO-C£[2-)  P^ak  integral  and  ^  is  as  defined  in  equation  (5). 
molar  ratio  TEOS  :  EtOH :  H2O :  HCl  =1:4;  0.8  :  0.0008 

6.  Mobility  of  silicates  in  TEOS  hydrolysis  and  condensation  as  determined  by  ^H  POSE 
NMR 

molar  ratio  TEOS  :  EtOH  :  H2O :  HQ  =  1 : 4  :  0.8  :  0.0008 
The  deviation  of  the  curve  is  estimated  to  be  ±  0.03  cm^/s, 
based  on  the  data  for  t  •>  180  h. 

*  molar  ratio  TEOS  :  EtOH :  H2O  =  1 : 4 : 0.8 

No  reaction  could  have  occurr^  since  diffusion  measurement  was  taken  within  0.5- 
1.0  hour  after  mixing. 

7.  ^^Si  NMR  spectra  ( ref.  TMS )  for  TEOS  reaction  mixture. 

molar  ratio  for  this  and  the  following  systems,  unless  stated  otherwise,  is 
TEOS  :  EtOH :  H2O :  HQ  =  1  : 4 : 0.8 : 0.0008 
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8.  Peak  assignments  of  NMR  multiplets  due  to  different  environment  of  adjacent 
silicon.  Symbols  follow  those  of  Kelts  and  Armstrong  [  24  ],  where  square  and  circle 
represents  resonating  and  adjacent  silicon,  respectively. 

9.  Time  evolution  of  silicon  species  as  detected  by  ^^Si  NMR. 

10.  Gas  chromatograms  for  TEOS  reaction  mixtures. 

1 1 .  Molecular  size  distribution  of  silicon  oligomers  as  determined  by  GC-FID. 

12.  Hydrolysis  and  condensation  of  TEOS:  water  and  ethanol  content,  calculated  from  ^^Si 
NMR  results. 


13.  ^^Si  NMR  spectrum  for  system: 

molar  ratio  TEOS  :  EtOH  :  H2O  :  HCl  =  1 : 4  :  0.5  :  0.0005 
reaction  time  1  month. 

1 4.  Mobility  of  silicon  oligomers  as  affected  by  the  extent  of  hydrolysis. 

Hydrolysis  initiated  by  addition  of  water  to  a  3-month-old  sample  with  the  same  molar 
composition  as  in  (13),  making  the  total  water-to-silicon  ratio,  W/Si,  to  be  0.8  and  1.0. 

1 5 .  Diffusion  coefficients  of  silicon  oligoimrs 
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